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COCHISE  Observations  of  Argon  Rydberg  Emission 

From  2  to  1 6  Micrometers 


1.  INTRODUCTION 
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Studies  of  the  radiative  and  collisional  properties  of  atomic  Rydberg  states 
have  generated  a  rich  literature  for  over  100  years.  1  A  great  many  of  the  more 
recent  experimental  investigations  have  employed  optical  excitation  methods  to 
form  highly  excited  states,  and  the  bulk  of  the  work  in  this  area  deals  with  alkali 
metals.  Relatively  little  attention  has  been  given  to  Rydberg  state  spectroscopy 
for  more  complex  atoms  such  as  O,  N,  and  Ar,  which  are  significant  constituents 
of  the  upper  atmosphere.  Furthermore,  in  high-altitude  electron  bombardments, 
as  in  low-pressure  laboratory  electric  discharges,  Rydberg  species  can  be  formed 
by  kinetic  processes  that  might  give  rise  to  non-Boltzmann  population  distributions 
among  the  upper  energy  levels.  Thus,  it  is  informative  to  study  the  excitation  of 
atmospheric  species  under  laboratory  conditions  representative  of  the  high  altitude 
environment.  We  report  here  the  first  direct  observations  of  long-wavelength 
infrared  (LWIR)  fluorescence  from  Rydberg  states  of  Ar  excited  in  a  low-pressure 
microwave  discharge,  as  studied  in  the  COCHISE  cryogenic  infrared  fluorescence 
facility. 

(Received  for  publication  4  August  1983) 

1.  Feneuille,  S. ,  and  Jacquinot,  P.  (1981)  Atomic  Rydberg  states,  Adv.  Atomic 
Molec.  Phys.  17:99. 
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Argon  is  a  significant  constituent  in  the  thermosphere,  having  a  vertical 

16  -2  2 

column  density  in  excess  of  10  cm  at  100  km.  In  a  strong  aurora,  there 
may  be  sufficient  excitation  of  high  Rydberg  levels  of  Ar  to  give  rise  to  measur¬ 
able  line  radiation  in  the  infrared,  where  radiative  lifetimes  of  such  transitions 
are  in  the  range  1  to  100  ps.  In  addition,  Ar  is  a  convenient  bath  gas  for  use  in 
all  COCHISE  infrared  chemiluminescence  studies,  and  its  spectral  signature  is 
present  to  some  extent  in  all  measurements.  For  these  reasons  we  have  con¬ 
ducted  a  detailed  investigation  of  the  Ar  emission  spectrum  observed  in  COCHISE. 

The  principal  objectives  of  this  investigation  are  to  characterize  the  Ar 
spectrum  from  2  to  16  pm  at  the  highest  possible  spectral  resolution,  and  to 
identify  as  fully  as  possible  the  sources  of  the  emission  and  the  kinetic  mech¬ 
anism  of  the  excitation.  Our  greatest  interest  lies  in  the  spectral  region  beyond 
8  pm,  where  the  data  are  unique  and  exhibit  an  unusual  spectral  distribution 
peaking  near  12  pm.  To  accomplish  these  goals,  we  (1)  obtained  considerable 
spectral  data  over  a  range  of  spectral  resolution  and  discharge  operating  condi¬ 
tions;  (2)  compared  the  observed  emissions  to  transitions  predicted  from  optical 
selection  rules;  and  (3)  compared  the  data  to  simulated  dipole-allowed  spectra 
computed  using  transition  probabilities  estimated  from  the  Coulomb  approxima¬ 
tion. 

For  Rydberg  atoms  in  which  a  single  electron  is  excited  to  an  orbital  of  large 
principle  quantum  number,  that  is,  near  the  point  of  ionization,  the  excited  elec¬ 
tron  moves  in  a  near-Coulombic  field  of  essentially  unit  charge.  Thus,  the 
transition  frequencies  in  Rydberg  spectra  are  often  represented  in  terms  of  a 
simple  hydrogenic  approximation.  In  COCHISE,  we  observe  infrared  transitions 
from  states  with  principal  quantum  numbers  4  <  n  <  12  (20,  000  to  1,  000  cm 
below  the  ionization  limit).  For  these  transitions,  we  find  that  the  hydrogenic 
approximation  does  not  apply,  even  if  quantum  defects  are  considered.  Further¬ 
more,  for  the  lower  values  of  n  exhibited  in  the  data,  whose  transitions  predom¬ 
inate  at  shorter  wavelengths  (X  <  5  ym),  we  observe  that  certain  optical  selection 
rules  break  down,  and  that  the  Coulomb  approximation  is  clearly  inappropriate 
for  predicting  spectral  intensities.  Finally,  the  spectral  distribution  at  long 
wavelengths  suggests  excitation  of  the  emitting  species  via  a  cascade  mechanism 
initiated  by  ionic  recombination. 


2.  Jacchia,  L.G.  (1977)  Thermospheric  Temperature.  Density,  and  Com¬ 
position:  New  Models,  Smithsonian  Astiophysical  Observatory  Special 
Report  375. 
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2.  EXPERIMENTAL  OBSERVATIONS 


The  design  and  operation  of  the  COCHISE  facility  have  been  described  in 
detail  elsewhere.  In  brief,  the  entire  radiative  environment  of  the  experi¬ 

ment  is  maintained  at  a  base  temperature  of  ~20  K,  which  effectively  eliminates 
background  radiation  within  the  2  to  20-^m  operating  range  of  the  apparatus.  The 
detection  system  consists  of  a  scanning  grating  monochromator  and  a  liquid- 
helium -cooled  Si:As  detector.  Reagent  gases  are  introduced  through  temperature- 
controlled  feed  lines  to  the  reaction  cell  shown  in  Figure  1.  A  flowing  Ar  mixture 
at  ~1  Torr  passes  through  four  parallel  microwave  discharges  (2450  MHz,  50  W) 

1.  1  cm  ID,  prior  to  expanding  into  a  low-pressure  (~3  mT),  cryogenieallv 
pumped  interaction  volume.  A  counterflow  gas  enters  the  volume  from  the  oppo¬ 
site  side  to  combine  with  the  discharged  gas  in  a  stagnation  region  near  the  axis 
of  the  cylindrical  field-of-view  of  the  detector.  In  some  cases  (for  example, 
Reference  8),  this  interaction  consists  of  a  chemical  reaction  under  nearly  single¬ 
collision  conditions.  In  other  experiments,  all  the  emission  observed  arises 

3  10 

from  processes  occurring  in  the  discharge  sidearms.  ’ 

In  typical  COCHISE  chemiluminescence  experiments,  Ar  acts  as  a  carrier 
gas  for  small  (<10  percent)  flows  of  reagents,  for  example,  Ng  or  Og,  which  are 
to  be  excited  in  the  microwave  discharges.  Such  conditions  typically  employ  Ar 


3.  Rawlins,  W.  T. ,  Caledonia,  G.E. ,  and  Kennealy,  J.  P.  (1981)  Observation 

of  spectrally  resolved  infrared  chemiluminescence  from  vibrationally 
excited  03<i/3),  J.  Geophys.  Res.  86:5247. 

4.  Rawlins,  W.T.,  Piper,  L.G.,  Green,  B.  D. ,  Wilemski,  G. ,  Goela,  J.  S. , 

and  Caledonia,  G.  E.  (1980)  LABCEDE  and  COCHISE  Analysis  II.  Vol.  I 
(Final  Report)  PSI-TR  207A,  Physical  Sciences  Inc.,  Andover,  Mass. 

5.  Rawlins,  W.  T. ,  Piper,  L.G.,  Caledonia,  G.  E. ,  and  Green,  B.  D.  (1981) 

COCHISE  Research,  PSI-TR  298,  Physical  Sciences  Inc. ,  Andover,  Mass. 

6.  Rawlins,  W.  T. ,  Murphy,  H.  C.,  Taylor,  R.  L. ,  and  Caledonia,  G.E.  (1980) 

Final  Report  on  COCHISE  Refurbishment.  PSI-TR  222,  Physical  Sciences 
Inc.,  Andover,  Mass. 

7.  Kennealy,  J.  P. ,  DelGreco,  F.  P. ,  Caledonia,  G.  E. ,  and  Rawlins,  W.T. 

(1979)  COCHISE:  Laboratory  spectroscopic  studies  of  atmospheric 
phenomena  with  high-sensitivity  cryogenic  instrumentation,  G.A.  Vanesse, 
Ed.,  Proc.  Soc.  Photo -Opt.  Instrum.  Eng.  191:151. 

8.  Kennealy,  J.  P. ,  DelGreco,  F.  P. ,  Caledonia,  G.  E. ,  and  Green,  B.  D.  (1978) 

Nitric  oxide  chemiexcitation  occurring  in  the  reaction  between  metastable 
nitrogen  atoms  and  oxygen  molecules,  J,  Chem.  Phys.  69:1574. 

9.  Caledonia,  G.  E. ,  Green,  B.  D. ,  Simons,  G.A.,  Kennealy,  J.P. ,  Robert, 

F.X.,  Corman,  A.,  and  DelGreco,  F.  P.  (1977)  COCHISE  Studies  I:  Fluid 
Dynamical  and  Infrared  Spectral  Analyses.  AFGL-TR-77-0281, 

AD  A053218. 

10.  Rawlins,  W.T.,  and  Armstrong,  R.A.  (1982)  COCHISE  Observations  of  O., 
Formed  by  Three-Body  Recombination  of  O  and  Q„.  AFGI.-TR-82-0273, 

AD  A123653. 
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Figure  1.  COCHISE  Reaction  Chamber.  The  cross-hatched 
region  delineates  the  field-of-view  of  the  detector  (not  shown) 


pressures  and  flow  rates  of  ~1  Torr  and  640  /nmol  s  \  respectively,  which  give 
residence  times  of  ~2  ms  in  the  chemically  active  portion  of  the  discharge.  The 
Ar  spectral  studies  reported  here  were  carried  out  for  pure  Ar  (Matheson,  UHP) 
at  discharge  pressures  of  0. 1  to  2  Torr  and  temperatures  of  80  to  150  K.  The 
approximate  discharge  residence  times  ranged  from  5  ms  at  0.  1  Torr  to  0.  8  ms 
at  2  Torr.  For  all  conditions,  a  large  number  of  sharp  spectral  features,  appar¬ 
ently  atomic  lines  and  multiplets,  were  observed  between  2  and  16  pm. 

The  observed  intensities  and  spectral  distributions  are  insensitive  to  the  flow 
rate  and  identity  of  the  counterflow  gas.  This  indicates  that  the  emission  arises 
not  from  the  field-of-view  but  from  the  discharge  tubes  themselves,  and  reaches 
the  field-of-view  by  multiple  reflections  within  the  reaction  chamber.  This  is 
consistent  with  our  expectation  that  the  emitting  species  can  be  excited  only  in 
the  active  discharge  and  that  their  radiative  lifetimes  (~1  to  10  ps  for  atoms) 
would  be  much  shorter  than  the  ~0.  5-ms  transit  time  from  the  end  of  the  discharge 
tubes  to  the  detector's  field-of-view.  We  expect  the  reduction  in  the  initial  line 
intensities  due  to  the  multiple  reflections  is  at  least  3  orders  of  magnitude  (that 
is,  of  order  10  reflections).  This  means  that  we  can  observe  only  the  most 
intense  features  of  the  many  that  must  arise  from  such  discharges,  that  is,  pri¬ 
marily  those  of  neutral  Ar. 


The  scattering  of  the  light  introduces  a  severe  restriction  into  the  interpre¬ 
tation  of  absolute  and  relative  intensities.  For  typical  chemiluminescence  exper¬ 
iments  in  COCHISE,  in-situ  blackbody  calibrations  permit  determinations  of  the 
absolute  and  relative  responsivity  of  the  optical  system  for  radiation  arising  from 
the  central  interaction  volume.  These  calibrations  have  been  applied  to  the  ob¬ 
served  Ar  spectra,  but  they  are  not  strictly  valid  for  scattered  light.  We  antici¬ 
pate  that  the  reflectivity  of  the  copper  wall  of  the  reaction  chamber  varies  rela¬ 
tively  weakly  and  uniformly  with  wavelength,  so  that  relative  intensities  over 
small  wavelength  intervals  (a  few  pm)  are  reasonably  representative  of  those 
emanating  from  the  discharge.  A  notable  exception  occurs  near  7  pm,  below 
which  the  sapphire  walls  of  the  discharge  tubes  will  transmit  the  plasma  radia¬ 
tion,  possibly  giving  artifically  large  intensities  relative  to  those  at  longer  wave¬ 
lengths.  We  present  our  observed  intensities  in  arbitrary  units  proportional  to 
W/cm  sr  pm  (as  determined  by  the  blackbody  calibrations)  to  give  a  meas¬ 
ure  of  the  apparent  relative  intensity  on  a  consistent  scale. 

Low -resolution  survey  spectra  of  the  observed  Ar-related  emission  are 
shown  in  Figures  2  to  4.  In  general,  there  are  many  intense  features  from 
2  to  5  pm,  a  striking  "window"  between  5  and  10  pm,  and  several  moderately 
strong  features  from  10  to  15  pm.  The  general  appearance  of  this  spectrum  is 
insensitive  to  pressure  and  temperature,  although  the  total  spi  ^tral  intensity 
diminishes  as  the  pressure  is  increased  (Tigure  5).  In  general,  as  N2  or  O 2  are 
admitted  to  the  discharge,  the  intensities  of  Ar-related  emissions  decrease. 

This  is  shown  for  the  case  of  02  in  Figure  6.  The  observed  signal  actually  passes 
through  a  maximum  near  1  percent  O 2,  suggesting  the  possibility  of  oxygen- 
related  emissions  within  the  ~0.  08-pm  (FWHM)  bandpass  of  the  measurement; 
unfortunately,  low  S/N  ratios  have  thus  far  prevented  high-resolution  scans  for 
these  discharge  conditions.  The  decrease  in  signal  upon  further  C>2  addition  is 
likely  due  to  quenching  of  the  emission  by  O  or  02-  The  pressure  effect  shown 
in  Figure  5  is  difficult  to  interpret  because  changing  the  Ar  pressure  in  the  dis¬ 
charge  sidearms  simultaneously  alters  the  flow  velocity  (that  is,  the  residence 
time)  and  the  E/N  in  the  discharge  plasma,  as  well  as  other  characteristics  of 
the  discharge. 

Most  of  the  detailed  spectral  studies  were  carried  out  at  0.  1  Torr,  where  the 
S/N  ratios  were  the  most  favorable.  Example  spectral  scans  at  medium  and  high 
spectral  resolution  are  shown  in  Figures  7  to  12.  All  the  features  exhibit  spectral 
characteristics  typical  of  atomic  lines  or  multiplets;  in  addition,  many  of  the  line 
groupings  cannot  be  separated  even  at  the  highest  resolution  possible  for  these 
measurements. 
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Figure  10.  Medium  Resolution  Spectrum  of  Ar  Emission 
10  to  12  nm,  0.  027-/jm  Resolution,  0. 1  Torr 
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Figure  11.  Medium  Resolution  Spectrum  of  Ar  Emission, 
12  to  15  Mm.  0.  027-jim  Resolution,  0.1  Torr 
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Figure  12.  High  Resolution  Spectrum  of  Ar  Emission  Near 
12.3  pm,  0.  013-pm  Resolution,  0.  1  Torr 


The  probable  sources  of  this  emission  are  Rydberg  states  of  Ar  or  Ar+ 


excited  in  the  discharge.  However,  the  levels  of  excitation  in  such  discharges 


are  rather  low  <[e -)  ~  1011  cm-3,  average  electron  energy  ~2  eV),  11  so  the 


total  Ar  number  density  is  several  orders  of  magnitude  smaller  than  that  of  the 
neutral  atom.  Thus,  the  observed  transitions  appear  to  arise  from  high  Rydberg- 
like  energy  levels  of  neutral  Ar.  However,  since  the  spectroscopic  literature 
does  not  extend  beyond  ~4  pm,  considerable  analysis  is  required  to  make  the 
spectral  assignments.  The  Rydberg  series  of  hydrogen-like,  n  —  n  -  1  transi¬ 
tions,  which  should  appear  in  the  spectra  of  all  highly  excited  atomic  gases,  is 
given  by 


0.09113 


n2(n  -  l)2 


2n  -  1 


where  \  is  the  line  position  in  pm  and  n  is  the  principle  quantum  number  of  the 
emitting  state  (n  2  5).  The  first  few  lines,  for  n  =  4,  5,  6,  7,  ...,  are  at  4.05, 
7.46,  12.36,  19.05,  ...  pm,  respectively.  However,  these  are  idealized  lines 
which  are  most  closely  approximated  by  transitions  involving  states  with  large 


11.  Kaufman,  F.  (1969)  The  production  of  atoms  and  simple  radicals  in  glow 
discharges,  Adv.  Chem.  Ser.  80:29. 
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orbital  angular  momentum  (high  i  -values);  other  transitions  at  nearby  wave¬ 
lengths  may  be  expected  for  Rydberg  states  with  lower  f  or  spin-orbit  coupling. 
The  Rydberg  wavelengths  of  4.  05  and  12.  36  pm  correlate  well  with  the  strong 
transitions  we  observe  at  ~4.  0  and  12.  3  pm;  however,  the  7.46-pm  transition  is 
conspicuously  absent  from  the  observed  spectra.  Furthermore,  the  observed 
spectra  are  clearly  complex  and  contain  many  strong  transitions  in  addition  to 
those  near  4  and  12  pm.  Indeed,  even  the  4-  and  12 -pm  features  are  extremely 
complex  when  examined  under  high  resolution  (refer  to  Figures  9  and  12).  Thus, 
in  order  to  interpret  the  spectral  data,  it  is  necessary  to  predict  all  possible 
transitions  according  to  the  energy  levels  and  optical  selection  rules  for  Ar,  and 
to  estimate  transition  probabilities  in  order  to  compare  the  experimental  spectra 
to  simulated  spectra  convolved  with  the  instrumental  resolution  function. 


3.  INFRARED  SPECTROSCOPY  OF  Arl 
3.1  Argon  Rydberg  States 

The  ground  state  configuration  of  the  Ar  atom  is  given  by 


ls22s22p63s23p6  . 


The  lowest-lying  excited  states  are  formed  by  promoting  a  p  electron  to  an  outer 
shell  to  give  configurations  of  the  form 


ls22s22p63s23p5  ni 


These  states  are  split  into  two  groups  of  levels  depending  on  the  state  of  the  core 

2 

electrons.  The  core  electrons  form  two  different  configurations,  P^2  and 
2p3/2*  w^ch  couple  t0  the  outer  electron  to  give  the  atomic  configuration. 

The  excited  states  are  described  using  the  intermediate  or  j-f  coupling 
scheme.  Specific  states  are  denoted  as 


nf[K]j  or  nl'[K]j  , 


where 


K  *  t  +  j  and  J  =  K  +  S 
Jc 


K  is  the  resultant  angular  momentum  formed  by  coupling  the  core  electron 

angular  momentum  j  to  the  outer  election  orbital  angular  momentum.  Primed 

states  denote  jc  »  1/2  and  unprimed  states  denote  Jc  *  3/2.  J  is  the  total  atomic 

momentum  formed  by  coupling  the  outer  electron  spin  argular  momentum,  S,  to 

K.  A  superscript  zero  is  used  to  denote  state  of  even  orbital  angular  momentum 

(f  ■  s,  d,  g).  The  primed  and  unprimed  manifolds  rapidly  divide  with  increasing 

n  into  two  distinct  groups  of  levels.  The  j  •  1/2  manifold  approaches  in  its 
o  c  2 

ionization  limit  the  P^2  states  Ar  ion;  jc  *  3/2  approaches  the  Pg^2 

state  ion.  The  respective  ionization  limits  of  the  j  =  3/2  and  j  =  1/2  manifolds 

-1  -i  c  c 

are  127109.9  cm  and  128541.3  cm  .  Energy  level  tables  have  been  compiled 

12  13 

by  Moore  and  more  recently  by  Bashkin  and  Stoner. 

A  simplified  energy  diagram  for  lowest  Ar  excited  states  is  shown  in 

Figure  13.  The  states  are  labeled  by  principal  quantum  number,  n.  The  energy 

span  for  states  of  given  n  and  I  is  depicted  by  a  vertical  line.  States  belonging 
o 

to  the  P jy g  core  are  marked  with  a  prime.  Two  sets  of  states.  4s  and  4s',  are 
below  the  scale  of  the  figure  and  therefore  are  not  shown.  One  immediate  obser¬ 
vation  is  the  rather  drastic  compression  of  states  for  higher  principal  quantum 

13 

numbers.  A  more  detailed  level  diagram  is  given  by  Bashkin  and  Stoner. 

3.2  Calculation  of  Radiative  Transition  Probabilities 

We  shall  summarize  the  calculation  of  radiative  transition  probabilities 
within  the  framework  of  a  one-electron  dipole  approximation.  The  selection 
rules  that  result  for  the  transition  between  states  njf  JKj  ]Jl  and  n2i2*K2*  J2 


1 1  "  l2  =  ±l 


K1  -  K2  =  0  (Kj  /  0).  ±1 


Jt  -  J2  *  0  (Jt  f  0).  ±1 


jc  ‘  jc  =  0 
C1  c2 


12,  Moore,  C.E.  (1971)  Atomic  Energy  Levels,  Vol.  1,  NSRDS-NBS35. 

13.  Bashkin,  S. ,  and  Stoner,  J.O.,  Jr.  (1978)  Atomic  Energy  -  Level  and 

Grotrian  Diagrams  2,  North-Holland  Publishing  Co.  ,  New  York. 


The  computed  line  strength  is  given  by 
S  =  -g-(2J1  +  1)(2J2  +  1)(2K1  +  1)(2K2  +  1) 


(  J1  J2 

|K2K1 


1 

1/2 


)  ( l2  K2 


<t 


2 


(2) 


where  the  bracketed  quantities  are  6-j  symbols  and  <i  Jj  T  || i  2 >  is  a  reduced 
matrix  element  of  the  dipole  operator.  It  is  expressed  as  the  following: 


IciJ  T||/2>|2  =  -^  (2^  +  l)(2i  2  +  l)C(f  ;  000)2 


00 

X  I  f  f.  Mi,  (r)rdr  | 2  ,  (3) 

J  1 1  l2 

o 

where  C(f  l,i2  ;  000)  is  a  Clebsch-Gordan  coefficient  and  it  (r)  and  f^  (r)  are 
radial  wave  functions  for  the  outer  electron.  We  will  use  the  Coulomb  approxima 
tion  for  the  evaluation  of  the  radial  integral. 

In  the  Coulomb  approximation  the  radial  wave  functions  are  approximated  by 
hydrogenic  wave  functions,  that  is, 

00  00 

I  s  J  f|  (r)rdr  =  J  ^°ui(r)  ^oul(r)rdr  ,  (4) 

12  12 
o  o 

where  is  a  hydrogenic  wave  function.  Several  methods  of  evaluation  of  I 

are  available. 

Bates  and  Damgaard14  and  others15'  16  have  presented  an  analytic  expres¬ 
sion  for  I: 


14.  Bates,  D.  R. ,  and  Damgaard,  A.  <1949)  The  calculation  of  the  absolute 

strengths  of  spectral  lines,  Phil.  Trans.  Roy.  Soc.  A242:101. 

15.  Statz,  H. ,  Tang,  C.  L. ,  and  Koster,  G.  F.  (1963)  Approximate  electromag¬ 

netic  transition  probabilities  and  relative  electron  excitation  cross 
sections  for  rare-gas  masers,  J.  Appl.  Phys.  34:2625. 

16.  Koster,  G.F. ,  and  Statz,  H.  (1964)  Probabilities  for  the  neon  laser  transi¬ 

tions,  J.  Appl.  Phys.  32 :2654. 
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1  c 

1 - TTTl - Ju+I  tr(nl+<  1+l)r(n1-i  1)r(n2+i2+l)r(n2-i2]_l/2 

(nj)  1  (i^)  * 


x  I  X 

kj-0  ^=0 


‘1  S  L(nl+n2>J  (  1 


+n2+2-k1-k2) 


The  principal  quantum  numbers  n^  and  n2  are  non-integers  and  are  evaluated  for 
each  state  from  the  following  relation 


n-  ( 


109737 


where  E..  is  the  ionization  limit  (cm  ;  of  the  state  manifold  and  E  is  the  term 
value  (cm  )  of  the  state.  The  coefficients  ak  and  a^  are  evaluated  from  the 
recurrence  relation 


na.  .  , 

ak  *  “Tfc  *  (n-k)(n-k+l)J  .  (7) 

with  aQ  *  1  for  each  set  of  coefficients.  The  summation  is  truncated  by  the  con¬ 
dition  kj  +  kj  <nj  +  n2  -1.  This  procedure  is  computationally  very  rapid  and 
works  well  for  nj  +  n2  <  using  double  precision  arithmetic.  The  problem  at 
higher  principal  effective  quantum  numbers  arises  from  cancellation  between 

large  numbers  in  the  summation  to  give  smaller  final  results.  For  higher  effec- 

17 

tive  principal  quantum  numbers  we  use  the  method  of  Katterbach. 

Katterbach  uses  a  semi-numerical  procedure  in  which  I  is  approximated  as: 


coul(r|,coul,rlrdr  „  y  ,coul, r„)*foul(r„l 

1  1  *0  M.  +  v  X.  n  V 

1  i  t^  =  0  1  i 


17.  Friedrich,  H. ,  Katterbach,  K. ,  and  Trefftz,  E.  (197  0)  On  the  accuracy 
of  machine  programs  for  calculating  o^  -'Hator  strengths  by  Coulomb 
approximation,  J.  Quant.  Spectrosc  ^  \  Trans.  10:11. 


The  hydrogenic  wave  functions  are  approximated  as 


^COul(r) 


[nT(n+l+l)r(n-i)) 


max 

wlf)  I  V- 


with  t  determined  by  the  inequality 
max 


t  <ntl<t  +1 
max  max 


and  the  at  defined  as  in  the  Bates -Damgaard  method.  The  step  size  Ar  is  gen¬ 
erally  taken  to  be  0.  125.  Larger  values,  Ar  =  0.  30,  give  identical  results  for 
the  cases  studied.  The  starting  point  rQ  is  taken  to  be  0.  25.  The  value  of  vmax 
is  determined  by  the  relation 


U.co(r„  >1 


max  |  4 


c°ul(r)| 


for  the  wavefunction  with  the  smaller  principal  quantum  number. 

We  have  found  that  these  two  methods  give  virtually  identical  results  in 
cases  when  comparison  is  possible.  The  Bates -Damgaard  method  is  more  com¬ 
putationally  efficient  and  therefore  used  when  applicable.  Other  cases  were 
computed  using  Katterbach's  method. 

The  Einstein  coefficient.  A,  of  the  emitting  species  is  inversely  proportional 
to  its  radiative  lifetime.  These  coefficients  are  related  to  the  wavelength- 
independent  atomic  linestrengths,  S,  by  the  relationship 

gA  .  l.my  io-*  yS  _  <„> 

XJ 

where  X  is  the  transition  wavelength  given  in  centimeters  and  g  the  degeneracy 
of  the  radiating  state. 

We  have  computed  gA  factors  for  radiating  states  from  n  =  4  to  n  =  22  for 
transitions  between  2  and  16  pm.  Additionally,  simulations  of  radiation  from 
hot  Ar  gas  have  been  computed.  The  results  of  these  calculations  are  presented 
in  the  next  section. 


4.  RESULTS  OF  SPECTRAL  ANALYSIS 


In  Table  1  we  list  tentatively  identified  Rydberg  transitions.  For  most 
observed  spectral  features  more  than  one  transition  is  given.  This  is  a  conse¬ 
quence  of  the  large  density  of  spectral  lines  for  excited  Ar  and  limited  spectral 
resolution  to  discriminate  between  the  various  transitions.  The  possible  transi¬ 
tions  were  determined  using  a  weaker  set  of  selection  rules  than  those  previously 
discussed  for  spectral  computations.  For  the  weaker  selection  rules  the  restric¬ 
tion  of  core  angular  momentum  invariance  is  dropped,  hence  transitions  between 

2  2 

the  Pj/2  anc*  ^3/2  manifolds  are  allowed.  Experimentally  these  transitions 
are  found  to  be  important  at  shorter  wavelengths,  for  example,  Humphreys18 
notes  many  core  changing  transitions,  several  quite  intense,  in  his  compilation 
of  observed  Ar  lines  between  1.2  and  4.  0  Mm. 

Table  1  includes  the  observed  peak  wavelength,  identified  possible  transi¬ 
tions,  and  predicted  wavelength  based  on  the  state  energies  listed  in  Moore's12 

1 3 

and  Bashkin  and  Stoner's  tables.  An  "H"  after  the  listing  indicates  that  the 
transition  is  listed  in  Humphreys'  compilation.  Table  1  includes  the  entire 
range  of  the  present  spectral  measurements,  from  approximately  2.  0  to  16.  0  Mm. 

Some  general  conclusions  about  the  types  of  identified  transitions  may  be 
made.  There  is  a  general  tendency  to  see  radiation  from  higher  energy  states, 
that  is,  higher  principal  quantum  numbers,  at  longer  wavelengths.  This  seems 
reasonable  since  the  energy  gaps  are  too  large  for  low  quantum  numbers  to 
radiate  in  the  LW1R.  However,  with  the  exception  of  the  peaks  observed  around 
12.2  to  12.4  Mm  (and  to  some  extent  around  2.4  Mm)  it  is  rarely  necessary  to 
invoke  states  of  quantum  numbers  greater  than  eight.  Two  possible  explanations 
are  that  excited  state  Ar  formation  does  not  produce  such  states  or  that  such 
states,  which  have  larger  radiative  lifetimes,  are  quenched  by  gas  phase  or 
wall  collisions. 

A  second  point  is  that,  as  mentioned  above,  core  changing  transitions  are 
quite  important  at  shorter  wavelengths.  Several  peaks  seem  to  be  directly 
attributable  to  such  transitions:  2.100,  2.212,  2.445,  2.634,  and  4.  116  Mm. 

The  significance  of  such  transitions  makes  the  use  of  the  Coulomb  approximation 
rather  dubious  since  they  are  beyond  the  scope  of  that  treatment.  It  seems  that 
for  longer  wavelengths  (8  to  16  Mm)  core  changing  transitions  play  a  less  signifi¬ 
cant  role.  The  Coulomb  approximation  may,  therefore,  be  more  reliable  for 
longer  wavelength  radiation.  (The  12.2  to  12.4  Mm  feature  is  an  exception.) 

This  conclusion  is  also  physically  appealing  since  for  higher  quantum  numbers 
the  outer  electrons  become  more  "hydrogen -like". 

18.  Humphreys,  C.  J.  (1974)  First  spectra  of  neon,  argon  and  xenon  136  in  the 
1.  2 -4.  0  Mm  region,  J.  Phys.  Chem.  Ref.  Data  2:519. 


24 


Table  1.  Tentatively  Identified  Spectral  Transitions  for  Argon 
Between  1.  95  and  16.  0  pm 


Table  1.  Tentatively  Identified  Spectral  Transitions  for  Argon 
Between  1. 95  and  16.  0  (Contd) 


Table  1.  Tentatively  Identified  Spectral  Transitions  for  Argon 
Between  1.  95  and  16.  0  pm  (Contd) 


Observed  Peak 

Predicted 

Wavelength 

Wavelength 

(pm) 

Transition 

(pm) 

3.  0236 
3.  0236 
3.  0223 
3.  1953 
3.  1995 
3.  2048 
3.2166 
3.  5068 
3.  5349 
3.5302 
3.  5280 
3.  5909 
3. 6971 
3. 6969 
3. 6920 
3.  9180 
3. 9330 
3. 9296 
3. 9363 
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Table  1.  Tentatively  Identified  Spectral  Transitions  for  Argon 
Between  1.  95  and  16.  0  um  (Contd) 


Observed  Peak 

Predicted 

Wavelength 

Wavelength 

(pm) 

Transition 

(pm) 

9p[l/2l1  -  5dt3/2l1 

3.  960 

9pll/2l2  -  5d[ 3/2 1  j 

3.  928 

6s[3/2]°  -  5p[3/2]2 

3. 9972 

6d'[3/2]J  -6p'[l/2]0 

4. 0069 

6s[3/2l2  -  5pl3/2)2 

4.  0001 

7p'[3/2l2  -  5dl7/2]° 

4.  0027 

6d'[3/2l2  -  6pl3/2l2 

4. 1129 

6d[3/2]°  -  6pl  1/2 ]Q 

4.2766 

5d[3/2]°  -4fl3/2)2 

4.2999 

5d(3/2)°  -4fl3/2)1 

4.2993 

7p[l/2]0  -*  6s'[1/2)° 

4. 3698 

5d'[5/2]g  -  6pl3/2]1 

4. 3730 

6s'(l/2l1  -  5p'[ l/2)0 

4. 3661 

8pl  1/2],  -  5dt 3 /2 lo 

4. 3680 

Table  1.  Tentatively  Identified  Spectral  Transitions  for  Argon 
Between  1.  95  and  16.  0  pm  (Contd) 


Observed  Peak 
Wavelength 
(pm) 

Transition 

Predicted 

Wavelength 

(pm) 

4.378 

4d'l5/2}°  -  5p[3/2)2 

4. 3775 

4d[3/2]°  -  5p[l/2)0 

4. 3768 

5d[3/2]°  -  4ft 5/2 12 

4. 3779 

6p[5/2l2  -  4dl5/2]° 

4. 3759 

4.440 

5d'(5/2]°  -  6p[3/2l1 

4.4478 

6dt  1/2]q  -  6p[3/2)1 

4.4410 

4.536 

5f'(5/2]  x  -*  5dll/2]° 

4. 5361 

4.720 

5p[ 5/2 J3  -  3d’[5/2]° 

4.7151 

5ft7/2l3  -  4d'[3/2l2 

4.7238 

6ft5/2]2  -  5dtl/2]J 

4.7198 

6d[7/2]g  —  4f'[5/2]2 

4.7185 

6d[7/2]°  -»  4f'[ 5/2 ]g 

4.7179 

6d[7/2]°  -4f'l7/2]3<4 

4.7158 

8d[5/2l2  —  7  p[  1  /  2  ]  1 

4.7184 

6.617 

6p[  3 / 2]  2  -  6st3/2]° 

6.  6205 

6p *[ 3/2 ]  j  -  6s[  l/2lg 

6. 6090 

7d'[3/2]°  -  7p[  1/2]  j 

6.6188 

10. 050 

4f[9/2]4  -  4d[7 / 2 ]g 

10.  052 

5d[ 7 /2  )g  -  6p[ 5/2 ]g 

10.  052 

10.330 

5dl7/2]°  -  6pl5/2l2 

10.  327 

-  V 
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Table  1.  Tentatively  Identified  Spectral  Transitions  for  Argon 
Between  1.  95  and  16.  0  (Contd) 


Observed  Peak 
Wavelength 
Oim) 

Transition 

Predicted 

Wavelength 

(uni) 

10. 840 

5d'[5/2]g  -  6p't3/2l2 

10.842 

11. 115 

4f'[7/2]3  4  -  4d't5/2]g 

11. 113 

4f't5/2]g  -  4d'[5/2]° 

11.  0889 

11.480 

5d'[5/2lJ  -  6p*[3/2]2 

11.498 

5d[7/2]°  -  6p[  5/2 lg 

11.484 

12. 010 

6dIl/2]°  -  6p'[3/2]2 

12.  006 

12. 130 

8dIl/2]°  -  8p[l/2)1 

12. 133 

12.  182 

13£13/2]1  2  -*9dl3/2]° 

12. 182 

18s[3/2}J  -  llp[3/2]2 

12. 186 

9s'[l/2]J  -8fl5/2]2 

12. 189 

12.229 

lld[3/2]°  -  8ft3/2l1  2 

12.230 

12.255 

5dl3/2]°  -  6pl3/2l2 

12.251 

12.318 

9dl5/2]°  -*7ft7/2]g  4 

12.319 

12.352 

12d[l/2]J  -  6f'[5/2]g 

12.355 

12p[5/2l2  g  -  8d[7/2lg 

12.348 

12.378 

10p[5/2]2  -  7d[5/2]° 

12.376 

12.385 

12d[5/2]2  -  10p[ l/2l1 

12.382 

13ft5/2l2  g  -  9d[7/2]g 

12. 388 

16dll/2]°  -  1  lp[  1  / 2 ] ^ 

12.389 

12.427 

10p(5/2]g  -7dt5/2)® 

12.426 

. * « •- 


Table  1.  Tentatively  Identified  Spectral  Transitions  for  Argon 
Between  1.  95  and  16.  0  pm  (Contd) 


Observed  Peak 
Wavelength 
(pm) 

Transition 

Predicted 

Wavelength 

(pm) 

13.060 

5f'[5/2l3  5d'[3/2l2 

13.  084 

5f't5/2l2  -  5d'[3/2]° 

13.  081 

7p[ 5/2)g  -  7s[3/2]° 

13. 060 

13.271 

5fl 3/ 2 ] x  -  5dll/2]° 

13.275 

5f[3/2]2  -  5d[l/2]J 

13.272 

13. 500 

5d[l/2]°  -  6p[5/2)2 

13.495 

7p(3/2l1  -  5dl3/2)° 

13.500 

13.665 

7p[l/2l1  -*  7 st 3/2 

13. 662 

14.625 

4fl7/2l3>4  —  4d[5/2l3 

14.619 

8s '[ l/2 ]q  -7p'tl/2]1 

14.626 

15.  570 

8p'[3/2J2  -  6d*[5/2]° 

15.578 

8s(3/2]°  -7p[3/2)2 

15. 584 

7d[5/2]g  -  6fl7/2J3 

15.611 

The  gA  factors  for  all  selection  rule  allowed  transitions  in  the  X  =  2-  to 
16-pm  wavelength  range  have  been  computed  for  initial  states  of  n  =  4,  5,  6,  7, 
8,  10,  12,  14,  16,  18,  20,  and  22.  The  results  of  these  calculations  are  pre¬ 
sented  in  Figures  14  to  25,  where  for  each  value  of  n,  gA  factors  are  plotted  vs 
transition  wavelength.  The  gA  values  are  normalized  to  the  largest  values  (for 
each  value  of  n)  within  the  2  to  16  pm  range.  The  maximum  value  of  gA  is  noted 
on  each  figure.  Also  noted  on  the  right-hand  vertical  axis  are  values  of  (gA)~*, 
proportional  to  the  predicted  lifetime  of  a  given  transition. 

Several  general  trends  in  the  values  of  gA  are  readily  seen.  Values  of  gA 
tend  to  decrease  as  the  transition  wavelength  increases.  This  behavior  is 
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dominated  by  the  «•  factor  in  Eq.  (11).  Values  of  gA  decrease  with  increasing 
values  of  principal  quantum  number,  n.  An  exception  to  this  rule  is  the  behavior 
of  n  =  4  and  n  =  5,  where  the  value  of  (gA)max  for  n  =  5  exceeds  that  for  n  =  4. 

The  largest  values  of  gA,  (gA)max>  vary  over  more  than  3  orders  of  magni¬ 
tude.  For  values  of  n  less  than  or  equal  to  6  there  are  many  lines  with  sub¬ 
microsecond  values  of  (gA)  However,  for  n  2  10  the  strongest  lines  are  in 
the  several  microsecond  lifetime  regime.  Note  that  these  time  scales  apply  to 
the  shorter  wavelengths,  2  to  4  /urn.  In  the  8  to  16  p m  range,  lifetimes  are  much 
longer.  In  that  range,  there  are  very  few  submicrosecond  values  of  (gA)"1  and 
only  for  n  ^  7.  For  n  >  10  there  are  no  allowed  transitions  with  lifetimes  less 
than  10  ps  and  for  n  >  14  (essentially)  none  with  lifetimes  less  than  100  ps. 
Consequently,  in  the  relatively  high  density  environment  of  the  COCHISE  appara¬ 
tus,  excited  states  with  high  quantum  values  of  n,  if  formed,  may  not  be  observed 
(between  8  and  16  pm)  due  to  more  rapid  gas  phase  or  wall  quenching  processes. 
This  point  will  be  discussed  in  a  later  section. 

These  calculations  were  made  within  the  framework  of  the  Coulomb  approx¬ 
imation,  usually  taken  to  be  order  of  magnitude  in  accuracy.  There  is  evidence 

that  we  may  have  some  confidence  in  the  general  accuracy  of  these  calculations 
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for  specific  transitions  in  Ar.  Recently,  Borge  and  Campos  '  have  measured 
lifetimes  for  several  ns[3/2]2  (n  =  6,  7,  8,  9)  states  and  compared  their  results 
to  computed  Coulomb  approximation  lifetimes  and  the  results  of  more  accurate 
calculations.  Their  comparisons  indicate  agreement  between  the  Coulomb  ap¬ 
proximation  predictions  and  measured  values  are  within  a  factor  of  2  for  these 
limited  transitions.  For  other  transitions,  we  have  found  the  agreement  to  be 
very  unsatisfactory. 

Calculations  were  performed  to  simulate  Ar  spectra  for  8  to  16  pm  as 
observed  in  COCHISE.  To  simulate  the  observed  spectra,  a  distribution  of 
populations  of  states  must  be  chosen.  We  have  chosen  a  Boltzmann  distribution 
with  a  temperature  of  10,  000  K.  Calculations  were  also  performed  at  tempera¬ 
tures  of  20,  000  and  50,  000  K.  They  showed  no  significant  variation,  since  the 
states  of  interest  are  at  high  energy  and  exhibit  only  slight  relative  population 
changes  with  temperature.  We  present  only  results  computed  at  10,000  K. 

The  use  of  a  Boltzmann  distribution  in  defining  state  population  is  convenient 
but  not  realistic.  In  the  disturbed  environment  of  a  microwave  discharge  it  is 
unlikely  that  equilibrium  is  achieved.  Specific  non-equilibrium  effects  have 


19.  Borge,  M.J.G.,  and  Campos,  J.  (1983)  Lifetimes  of  ns[3 /2  ]  2  (n  =  6, 7 ,  8,  9) 
levels  and  transition  probabilities  of  4p  -  ns  lines  of  Arl,  Phys.  Rev. 
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previously  been  observed  by  various  workers  in  Ar  discharges  and  should 
be  expected  here. 

In  Figure  26  are  the  predicted  spectra  from  excited  Ar  at  10,  000  K  and  the 
experimentally  determined  COCHISE  spectra  in  the  8  to  16  pm  range.  The  theo¬ 
retical  spectrum  was  computed  using  a  triangular  slit  function  of  width  0.  08  pm 
to  simulate  the  resolution  of  the  experimental  spectrum.  The  theoretical  inten¬ 
sity  is  plotted  in  relative  units  of  spectral  radiance.  Note  that  the  ordinate  scale 
is  linear.  The  energy  radiated  at  a  given  wavelength  is  given  by  the  formula 


-e./kT 

I(X)  =  £  giAi_je  a(X-X. .)e..  , 


where  g.A.  .  is  the  upper-state  degeneracy  times  the  Einstein  coefficient  for 

1  -Ei/Vr 

the  transition  between  states  i  and  j,  e  the  Boltzmann  factor  for  the  upper 

state,  6^  the  energy  of  the  transition,  and  ct(X-X_)  the  triangular  slit  function. 
The  sum  is  over  all  possible  transitions  allowed  in  the  dipole  approximation. 

The  comparison  between  theory  and  experiment  is  mixed.  In  the  interval 
8.  5  to  10  pm  very  little  radiation  is  measured  where  the  calculations  predict 
substantial  emission.  The  calculations  also  do  not  predict  the  shape  and  relative 
magnitude  of  the  large  spectral  feature  around  12.4  pm.  However,  with  these 
exceptions  the  theoretical  predictions  correlate  well  with  the  observed  peak 
positions  and,  within  factors  of  2  to  3,  with  relative  peak  heights  over  the  rest 
of  the  spectrum,  that  is,  the  theoretical  predictions  seem  to  agree  best  with 
experiment  in  the  10  to  12  pm  and  12.6  to  16.  0  pm  ranges.  There  are  several 
possible  explanations.  The  most  obvious  is  that  the  Coulomb  approximation  is 
failing  in  the  8.  5  to  10  pm  region.  It  is  curious,  however,  that  several  of  the 
observed  transitions  in  the  10  to  12  pm  or  12.  6  to  16.  0  pm  region  originate  from 
similar  upper  states  as  the  unobserved  transitions  in  the  8.  5  to  10  pm  region. 
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This  may  indicate  that  there  is  a  non -equilibrium  population  due  to  selective 
excitation  or  cascade  processes  in  the  radiating  gas. 

Identified  transitions  in  the  10  to  12  pm  and  12.  6  to  16.  0  pm  regions  are 
listed  in  Table  1.  Several  features  are  noteworthy.  None  of  the  identified 
transitions  involve  core  changing  transitions.  This  is  quite  unlike  the  identified 
transitions  at  shorter  wavelengths,  where  several  strong  core  changing  transi¬ 
tions  are  noted.  There  is  also  a  tendency  to  see  transitions  from  relatively  high 
excited  states.  Transitions  from  n  =  5  are  most  common  with  some  potential 
transitions  from  n  =  8  noted.  Transitions  from  higher  states  are  not  needed  to 
account  for  these  observed  spectra.  It  may  be  that  higher  states  are  not  pro¬ 
duced  or  that,  when  formed,  they  are  quenched.  A  discussion  of  these  possibili¬ 
ties  is  given  in  a  later  section  on  the  excitation  mechanisms  and  kinetics. 

The  strongest  predicted  transitions  in  the  8.  5  to  10.  0  pm  region  are  listed 
in  Table  2.  The  upper  states  of  these  transitions  are  reasonably  similar  to 
identified  transitions  noted  in  Table  1.  We  have  examined  this  group  to  see  if 
transitions  arising  from  the  same  upper  state  appear  in  other  portions  of  the 
spectrum.  Most  do  not  appear  as  singly  identified  transitions,  but  several 
appear  identified  with  other  transitions  that  may  well  cause  the  observed  feature. 
Only  transitions  from  the  4f ’ [7/2 ]  3 ,  4f *[5/2  J  3,  and  5d[3/2 ]  2  states  seem  to  be 
definitely  identified  at  11.115,  11.115,  and  12.255  pm,  respectively. 

For  the  strong  transitions  predicted  between  8.  5  to  10  pm,  we  might  specu¬ 
late  that,  for  states  not  radiating  elsewhere  in  the  observed  spectrum,  non- 
equilibrium  effects  from  selective  excitation  or  cascade  may  cause  them  not  to 
be  populated.  The  seemingly  missing  predicted  radiation  from  4f  states  between 
9.497  and  9.  671  pm  might  well  be  due  to  the  breakdown  of  the  Coulomb  approxi¬ 
mation  for  4f  states.  The  Coulomb  approximation  limits  are  generally  accepted 
to  be  n  -  f  2  1,  so  that  this  case  is  at  best  borderline.  The  possibility  of  non¬ 
equilibrium  populations  is  quite  intriguing  although  we  readily  admit  to  the  quan¬ 
titative  limitations  of  the  Coulomb  approximation. 

The  predicted  and  experimentally  measured  spectra  for  Ar  radiation  between 
11.  90  and  12.  50  pm  are  shown  in  Figure  27.  The  theoretical  spectrum  was  com¬ 
puted  using  a  triangular  slit  function  of  width  0.  013  pm  to  simulate  the  resolution 
of  the  experimental  spectrum.  The  theoretical  spectrum  is  again  plotted  in 
relative  units  of  power  radiated/wavelength.  Comparison  between  the  theoretical 
and  experimental  spectra  is  nominal.  Most  of  the  peaks  in  the  measured  spec¬ 
trum  are  accounted  for  except  for  those  around  12.35  to  12.1J  pm.  Identified 
lines  in  this  spectral  region  are  listed  in  Table  1.  Note  that  no  particularly 
strong  lines  are  predicted  at  12. 352  and  12.  385  pm.  The  origin  of  these  intense 
lines  is  not  apparent.  However,  the  observed  features  between  12  and  12.5  pm 
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Table  2.  Spectral  Transitions:  8.  50  to  10.  00  pm 


Transition 

Predicted 

Wavelength 

(pm) 

4f'l7/2]3  -4d'[5/2]J 

9.668 

4f'[5/2]3  ~4d'[3/2]° 

9.493 

5d£ 3/2 Jg  -  6p[3/2]1 

9.822 

5dl5/2]°  -*  6pl  1/2 lx 

9.757 

5d'[3/2]J  -  6p'tl/2]0 

9.7  57 

5d[3/2]J  -  6p[l/2l0 

9.  579 

5d[5/2]3  -  6p[3/2l2 

9.443 

5d[5/2]2  -  6p[5/2]2 

9.  172 

5d[5/2l3  -  6p[5/2)3 

8.  589 

6d'[5/2l3  -  5fl7/2l4 

9.  859 

6d[7/2l3  -  5f( 9/2l 4 

9.275 

6 dt 5/2] 2  -  5fl7/2]3 

9.  018 

7f'[5/2l2  -  Gd't3/2]J 

9.913 

7ft7/2]4  -  6d(5/2l3 

9.650 

4p[5/2l2  -  5d[7/2]° 

9.429 

7p'[5/2l1  -  5d'(5/2]2 

8.787 

7sl3/2l2  -  6p[3/2l2 

8.  551 

7p[5/2)3  -  5d[7/2]J 

8.  548 

can  be  accounted  for  by  invoking  transitions  from  extremely  high  levels,  n  >  9. 
The  difficulty  here  is  that,  because  of  the  relatively  long  radiative  lifetimes  for 
such  transitions,  large  steady-state  populations  of  these  levels  are  required  to 
explain  the  intense  emission  we  observe.  This  point  will  be  explored  further  in 
Section  5. 

In  summary,  the  theoretical  spectra  (computed  in  the  Coulomb  approxima¬ 
tion)  give  variable  agreement  with  the  LWIR  COCHISE  (8  to  16  pm)  spectra. 
Several  significant  discrepancies  exist,  and  specifically  the  theory  fails  to  pre¬ 
dict  the  relatively  strong  radiation  at  around  12.40  pm. 


5.  DISCUSSION 


5.1  Kinetic  Interpretation* 

Careful  inspection  of  the  transitions  listed  in  Table  1  leads  to  several 
interesting  observations.  For  convenience,  the  observed  manifold -to -manifold 
transitions  are  depicted  in  Figures  28  to  30. 

First,  emission  is  observed  from  states  up  to  within  a  few  hundred  cm-1  of 

the  ionization  limit  for  the  j  =  3/2  (unprimed)  core  configuration.  This  is  far 

higher  in  the  energy  manifold  then  has  been  probed  in  other  electron  bombard- 
2  1 

ment  experiments,  in  which  only  emission  in  the  visible  was  monitored. 

Second,  emitting  states  above  124,  000  cm  *  are  not  uniformly  represented 
by  observable  transitions.  This  is  contrary  to  what  we  might  expect  for  excita¬ 
tion  solely  by  direct  electron  impact.  Electron  impact  excitation  will  occur  in 
the  tail  of  the  electron  energy  distribution  with  cross  sections  similar  to  that  for 
ionization.  This  should  give  rise  to  a  uniform  population  distribution,  with  the 
relative  state  populations  decreasing  with  increasing  energy.  Collisional  relaxa¬ 
tion  may  modify  this  distribution  somewhat,  but  this  modification  should  not 
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strongly  affect  the  distributions  under  our  conditions.  Quenching  studies 

for  the  4p  and  4p'  manifolds  yield  collisional  deactivation  rate  coefficients  in  the 
_ i i  3-1 

10  cm  s  range,  probably  occurring  via  curve  crossings  on  the  Ar9  potential 

24  1  fi  « 

surface.  For  our  conditions,  0.  1  Torr,  100  K,  [Ar]  *  10  cm  ,  the  colli¬ 
sional  lifetime  of  Ar*  would  be  ~  10  ms,  that  is,  on  the  order  of  the  infrared 
radiative  lifetimes  for  the  higher  levels  (n  ~  10).  Additional  rapid  quenching  by 
electronic -to -translational  energy  transfer  may  be  expected  for  states  above 
n  ~  12,  where  the  energy  defects  are  on  the  order  of  kT  or  less.  Thus,  we 
expect  quenching  to  deplete  levels  above  125,  000  cm  1,  but  in  a  non-selective 
way.  The  implication  is,  then,  that  the  few  observed  emitting  states  above 
125,  000  cm  ,  which  radiate  predominantly  at  long  wavelengths,  are  directly 
populated  by  some  rapid  kinetic  process.  This  is  consistent  with  the  results  of 
the  spectral  analysis,  which  show  that  the  upper-level  electronic  population  dis¬ 
tribution  cannot  be  described  by  a  Boltzmann  expression,  but  rather  exhibits  a 
population  enhancement  for  LWIR -radiating  states.  The  LWIR  emissions  from 
these  states,  as  well  as  their  collisional  deactivation,  appear  to  be  the  initial 
steps  in  a  (largely  radiative)  cascade  sequence  that  populates  the  lower  levels. 

As  the  cascade  moves  down  the  energy  "ladder",  the  radiative  lifetimes  become 
shorter,  quenching  becomes  less  important,  and  the  radiative  transitions  occur 
at  shorter  wavelengths.  Below  —  122,000  cm  *,  virtually  every  electronic  state 
manifold  can  be  accounted  for  in  the  infrared  spectrum.  Note  that  the  relative 
populations  within  manifolds  can  still  be  non-Boltzmann  due  to  the  selectivity  of 
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the  radiative  transitions  in  and  out.  This  effect  has  been  observed  even  at  high 
22  23 

pressures.  *  Such  selective  excitation  within  manifolds  can  explain  the  ob¬ 
served  absence  of  certain  predicted  transitions,  as  discussed  in  the  previous 
section. 

Third,  the  observed  transitions  occur  primarily  in  the  j  =  3/2  manifold. 

Of  the  few  j  =  1/2  transitions  we  observe,  most  of  them  occur  at  shorter  wave¬ 
lengths,  and  the  highest  states  observed  lie  at  substantially  lower  energies  than 
those  observed  for  j  =  3/2.  The  implication  here  is  that  there  is  little  or  no 
direct  formation  of  Ar*  with  the  j  =  1/2  core,  and  that  the  higher  states  we  ob¬ 
serve  are  populated  initially  by  collisional  energy  transfer  from  Ar*  (j  =  3/2). 

Fourth,  there  are  some  manifolds  at  intermediate  points  in  the  cascade  for 
which  there  are  several  radiative  inputs  and  few  radiative  outputs.  This  is 
especially  true  for  the  6p  manifold,  and  also  to  some  extent  for  5d,  4f,  4d,  and 
5p.  The  "bottleneck  effect"  in  evidence  for  these  manifolds  may  be  expected  to 
give  unusual  state  distributions  within  each  manifold  and  intense  radiative  transi¬ 
tions  at  selected  wavelengths,  especially  in  the  near  infrared  (X  <  2  Aim),  where 
the  strongest  emissions  for  these  states  will  occur. 

It  remains  to  identify  the  excitation  process  responsible  for  initiating  the 
cascade.  An  obvious  candidate  is  the  dissociative  recombination  of  Ar2 

e"  +  ArJ  -»  Ar  +  Ar*  ,  (12) 

where  Ar*  is  in  one  of  many  possible  Rydberg  states.  Direct  recombination  of 
Ar+  with  electrons  in  a  three-body,  bielectronic  process  and  will  be  extremely 
slow  in  a  low-pressure,  weak  plasma,  such  as  occurs  in  our  discharges.  Simi¬ 
larly,  two-body  radiative  dielectronic  recombination  of  Ar+  will  not  be  significant 

under  our  conditions.  Reaction  (12)  has  been  carefully  studied  by  Shiu  and 
2 1 

Biondi  for  electron  temperatures  of  300  to  7500  K.  By  monitoring  radiation 
in  the  visible,  they  observed  formation  of  Ar  (4p,  4p')  in  the  recombination  of 
Ar2  with  300  K  electrons,  and  formation  of  states  as  high  as  5d  when  7200  K 
electrons  were  used.  They  deduced  that  the  observed  states  were  populated 
directly  from  Ar2,  whose  zero-point  energy  lies  near  117,000  cm  1  (binding  en¬ 
ergies  of  1.4,  1.23  ±  0.02,  and  1.33  ±  0.02  eV  are  reported  in  References  25  to  27, 

25.  Mullikan,  R.S.  (1970)  Potential  curves  of  diatomic  rare-gas  molecules  and 

their  ions,  with  particular  reference  to  Xe2*.  J.  Chem.  Phys.  52:5170. 

26.  Ng,  C.  Y. ,  Trevor,  D.  J. ,  Mahan,  B.  H. ,  and  Lee,  Y.  T.  (1977)  Photoioniza¬ 

tion  studies  of  the  Kr?  and  Ar0  van  der  Waals  molecules,  J.  Chem.  Phys. 
6i6:446.  ‘  "  . 

27.  Moseley,  J.  T. ,  Saxon,  R.P.,  Huber,  B.A.,  Cosby,  P.  C. ,  Abouaf,  R.,  and 

Tadjeddine,  M.  (1977)  Photofragment  spectroscopy  and  potential  curves  of 
ArJ,  J.  Chem.  Phys.  67 :1659. 
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respectively).  However,  since  infrared  emission  was  not  monitored  in  their 

experiments,  it  is  not  clear  whether  the  observed  states  were  populated  directly 

or  Dy  radiative  cascade  from  above.  In  any  case,  Shiu  and  Biondi's  experiments 

show  that  energetic  electrons  can  recombine  with  Ar^  to  form  Ar*  with  most  of 

the  excess  energy  going  into  electronic  excitation  of  the  Ar*.  For  the  microwave 

4 

discharge  case,  with  an  electron  temperature  of  2  -  3  X  10  K,  one  can  expect 
population  of  states  near  the  ionization  limit,  as  we  observe  spectroscopically. 
(There  may  also  be  substantial  formation  of  Ar  (  Pg/g)  by  this  route* )  Of  course, 
we  cannot  eliminate  the  possibility  that  intermediate  states  of  Ar*  are  directly 
populated  as  well,  which  will  exacerbate  the  non-Boltzmann  character  of  the 
cascade  process. 

Reaction  (12)  can  only  explain  our  observations  if  there  is  sufficient  forma¬ 
tion  of  Arg  in  a  microwave  discharge.  Arg  can  be  formed  by  three-body  recom¬ 
bination  of  Ar+: 


Ar+  +  Ar  +  Ar  —  Arg  +  Ar  .  (13) 

1  e  _ 1 

However,  Reaction  (13)  has  a  rate  coefficient  of  the  order  of  10  cm  s  at 
2  8 

room  temperature,  giving  it  an  effective  bimole cular  rate  coefficient  of 
-15  3-1 

•'•10  cm  s  at  0. 1  Torr.  By  comparison,  direct  electron  impact  excitation 

of  Ar*  should  occur  with  an  effective  rate  coefficient  similar  to  that  for  ioniza- 

29  -12±1 

tion;  Caledonia  has  calculated  the  latter  quantity  to  be  in  the  range  10 

3-1  + 

cm  s  for  typical  microwave  discharges.  Thus,  we  anticipate  that  the  Arg 

mechanism,  with  Reaction  (13)  as  a  rate-limiting  step,  would  not  be  able  to 

compete  with  electron  impact  excitation  as  a  source  for  Ar*  unless  thf  pressure 

2 1 

is  considerably  larger.  However,  as  pointed  out  by  Shiu  and  Biondi,  the 

reverse  of  Reaction  (12),  associative  ionization  of  Ar*.  is  probably  very  facile, 

-9  3-1 

with  a  rate  coefficient  of  at  least  10  cm  s  .  llius ,  it  is  possible  to  form 
Ar„  from  less  energetic  Ar*  produced  in  electron  impact,  and  then  recombine 
the  Ar2  with  "hot"  electrons  to  form  highly  energetic  Ar**: 


e  +  Ar  —  Ar*  +  e 


28.  Good,  A.  (1975)  Third-order  ion-molecule  clustering  reactions,  Chem. 

Rev.  75:561. 

29.  Caledonia,  G.E.  Unpublished  work.  See  Appendix  C  of  Rawlins,  W.  T. , 

Piper,  L.G.,  Gelb,  A.,  Lucht,  R.A.,  and  Caledonia,  G.E.  (1983) 
COCHISE  Research,  PSI-TR  349,  Physical  Sciences  Inc. ,  Andover,  Mass. 
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Ar **  +  Ar 


(16) 


ehot  +  Ar+2  ~ 


Here  denotes  the  relatively  energetic  "average"  electron  in  the  plasma. 
Although  the  recombination  coefficient  decreases  with  increasing  electron  energy 
the  electron  energy  distribution  is  such  that  most  recombinations  will  occur  with 


relatively  energetic  electrons.  Ar**  is  then  a  highly  excited  Rydberg  species 


with  a  P„ 


*3/2  core  configuration,  and  is  subject  to  collisional  and  radiative 

deactivation  as  outlined  above,  as  well  as  to  rapid  electron-impact  ionization. 

2 

The  specificity  of  the  process  for  the  P3/2  core  can  understood  in  terms  of 

the  adiabatic  dissociation  of  the  Art  —  e“  associated  pair  along  the  Ar  — 

+  2  *  ^  + 

Ar  (  P^/2^e  coordinate:  the  complex  behaves  like  Ar2  and  dissociates  to  form 

the  lowest-energy  ionic  (core)  configuration  of  the  atomic  "ion". 

Reactions  (15)  and  (16)  in  effect  set  up  a  quasi -equilibrium  between  Ar2  and 


the  atomic  Rydberg  manifold,  such  that  the  net  removal  of  Rydberg  atoms  is 
controlled  primarily  by  radiative  and  collisional  cascade.  Recombination  of 
Ar2  and  e”ojd  (the  reverse  of  (15))  can  also  occur;  however,  this  reaction  only 
feeds  back  Ar*  and  should  not  constitute  a  major  loss  term  for  Ar2.  Since  this 
reaction  is  dominated  by  ambient-temperature  electrons,  we  anticipate  that,  in 
the  active  discharge,  the  principal  fate  of  Ar2  is  recombination  with  e^ot  in 
Reaction  (16).  In  steady  state,  the  Ar2  and  Ar**  number  densities  are  given  by: 


[Ar+]  =  [Ar] 


kc[ehotJ 


(17) 


[Ar**] 


k  [e'][Ar] 


(18) 


where  L  is  the  effective  first-order  loss  rate  of  Ar**  by  collisional  and  radiative 

6*1 

deactivation  and  is  of  order  10  s  .  We  can  see  from  Eq.  (18)  that  [Ar**] 
depends  upon  the  effective  rate  for  electron -impact  excitation  of  Ar*.  This  rate 
is  extremely  sensitive  to  pressure  through  the  field/density  ratio  E/N  for  the 
discharge;  indeed,  the  ionization  rate  coefficient  increases  2  orders  of  magni- 
tude  for  a  two-fold  increase  of  E/N.  Thus,  the  pressure  effect  we  observe  is 
undoubtedly  due  largely  to  changes  in  discharge  conditions  such  as  E/N  and  [e  ]. 


10  -3 

An  order -of -magnitude  estimate  from  Eq.  (18)  gives  [Ar**]  ~  10  cm  in 

4  -1 

the  active  discharge.  For  a  gA  factor  of  ~  10  s  ,  this  value  yields  a  spectral 
-5  -2  -1  -I 

radiance  of  —  10  W  cm  sr  nm  at  12  pm.  This  is  to  be  compared  to  the 

-9  -2  -1  -1 

apparent  value  of  ~  2  x  10  W  cm  sr  pm  determined  in  our  experiments 
by  neglecting  geometric  effects.  Given  the  three -order -of -magnitude  reduction 
in  signal  we  expect  due  to  multiple  light  scattering,  the  predicted  and  observed 
emission  levels  are  consistent  with  each  other.  Clearly,  direct  observation  of 
the  discharge  plasma  in  COCHISE  would  provide  a  great  deal  more  information; 
however,  such  experiments  require  extensive  reconfigurations,  which  are  pres¬ 
ently  incompatible  with  our  on-going  chemiluminescence  investigations.  We  are 
considering  future  experiments  incorporating  such  modifications. 

5.2  Upper  Atmospheric  Implications 

Neutral  Ar  is  a  significant  constituent  of  the  thermosphere,  having  mixing 
-2  -3  2 

ratios  of  10  at  90  km  and  10  at  140  km.  The  number  density  of  Ar  is  con¬ 
sistently  about  3  or  4  orders  of  magnitude  larger  than  that  of  0„  in  this  altitude 

_3 

range.  At  the  low  pressures  characteristic  of  these  altitudes  (<  10  Torr), 
excitation  of  Ar  Rydberg  emission  in  an  electron  event  could  occur  via  the  same 
general  mechanism  that  appears  to  control  its  production  in  the  COCHISE  experi¬ 
ments,  with  an  important  exception:  deactivation  of  Ar*  and  Ar**  would  be 
entirely  radiative.  This  means  that  much  of  the  Ar*  formed  initially  by  electron- 
impact  would  radiate  rather  than  associatively  ionize  with  Ar.  Thus,  we  might 
expect  the  radiation  signature  to  be  more  characteristic  of  direct  excitation  by 

electrons,  that  is,  the  LWIR  component  would  be  diminished.  In  a  full  limb- 

4  -1 

viewing  aspect,  for  gA  —  10  s  near  10-12  pm,  a  typical  instrumental  noise 

level  of  ~  10  W  cm  ^  sr  1  would  correspond  to  a  minimum  detectable  emitter 

-3  4  -3  -1 

concentration  of  order  1  cm  ,  or  a  production  rate  of  Ar*  of  ~  10  cm  s 

30 

Judging  from  predicted  auroral  ion-pair  production  rates,  this  level  of  direct 


Ar*  formation  is  not  unreasonable  for  a  strong  (IBC  III)  aurora.  However,  in  a 
normal  atmospheric  scenario,  it  is  unlikely  that  sufficient  excitation  could  occur 
to  give  measureable  emission  for  zenith  aspects,  which  require  number  densities 
2  orders  of  magnitude  larger.  For  extraordinary  events  such  as  nuclear  phe¬ 
nomena,  Ar  Rydberg  emission  near  12  pm  (and  elsewhere  in  the  LWIR)  would 
probably  be  significant  for  any  viewing  geometry. 

The  COCHISE  results  for  Arl  may  hold  implications  for  atmospheric  OI  and 

+  +  + 

NI  emission  in  the  LWIR.  Unlike  Ar2,  02  and  N2  are  readily  formed  by  direct 
ionization  of  air;  the  recombination  of  these  ions  with  energetic  electrons  may 

30.  Caledonia,  G.  E.  ,  and  Kennealy,  J.P.  (1982)  NO  infrared  radiation  in  the 
upper  atmosphere.  Planet.  Space  Sci.  30:1043. 
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lead  to  the  formation  of  high  Rydberg  states  of  O  and  N.  If  the  initially  populated 
Rydberg  states  are  high  enough,  substantial  line  radiation  could  occur  in  the 
10-  to  16-jmi  region  in  a  strong  aurora.  However,  the  spectral  signature  is 
difficult  to  predict  theoretically  because  of  the  non-Boltzmann  character  of  the 
initial  formation  and  cascade.  We  plan  to  pursue  this  point  with  laboratory 
experiments  on  low-pressure  (<  1  Torr)  Ar/Ng  and  Ar/Og  discharges  and  with 
spectral  simulations  using  the  methodologies  we  have  applied  to  the  Arl  data. 

6.  CONCLUSIONS 

We  have  used  the  COCHISE  apparatus  to  obtain  a  new  and  unique  set  of 
spectroscopic  data  on  infrared  emission  from  Ar  Rydberg  atoms  by  observing 
scattered  light  from  a  low-pressure  microwave  discharge  plasma.  These  meas¬ 
urements  are  the  first  laboratory  observations  of  LWIR  Rydberg  emission,  and 
exhibit  a  rich  spectrum  from  2  to  16  ^m.  Commonly  used  hydrogenic  approxima¬ 
tions  often  applied  to  other  Rydberg  atoms  cannot  account  for  the  complex  nature 
of  this  spectrum,  but  all  the  observed  transitions  can  be  identified  as  Arl  lines 
or  groupings  using  tabulated  energy  levels  together  with  the  strong  (dipole- 
allowed)  selection  rules  and,  in  several  cases,  also  allowing  the  core  configura¬ 
tion  to  change  in  the  transition.  A  Coulomb  approximation  to  the  dipole  moment 
function  is  only  partially  successful  in  accounting  for  the  observed  intensities  of 
the  dipole -allowed  transitions,  as  it  fails  to  predict  detailed  spectral  distribu¬ 
tions  at  shorter  wavelengths  (A  <  5  fjm),  where  core-changing  transitions  are 
important.  The  spectral  distribution  in  the  LWIR  portion  of  the  spectrum  exhib¬ 
its  a  remarkable  peak  near  12  um  that  can  only  be  accounted  for  by  a  non- 
Boltzmann  population  distribution  with  a  strong  enhancement  of  very  high  states 
within  a  few  hundred  cm  *  of  the  ionization  limit.  The  most  likely  source  of 
this  distribution  is  the  dissociative  recombination  of  Ar*  with  energetic  electrons 
in  the  active  discharge,  followed  by  rapid  collisional  and  radiative  cascade. 

These  infrared  measurements  are  the  first  observations  of  such  high  excitation 
in  dissociative  recombination;  previous  workers  have  studied  only  visible  emis¬ 
sions,  which  are  the  end  results  of  the  cascade  process. 

The  Arl  emission  studied  here  may  have  some  bearing  on  high-altitude  LWIR 
observations,  particularly  near  12  um.  However,  although  Ar  is  a  significant 
thermospheric  species,  it  is  difficult  to  form  substantial  levels  of  Ar2  at  high  al¬ 
titudes,  so  Arl  radiation  is  probably  observable  only  in  limb-viewing  aspects  of 
relatively  strong  electron  disturbances.  Of  greater  atmospheric  significance  are 
the  implications  of  these  results  for  ionized  air,  that  is,  that  Og  and  N2  may  also 
recombine  with  energetic  electrons  to  initiate  an  infrared  cascade  sequence  of 
OI  and  NI  Rydberg  emission.  We  will  pursue  this  latter  issue  in  further  labora¬ 
tory  experiments  on  the  COCHISE  facility. 
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